Practicality of spin chain 'wiring' in diamond quantum technologies 



(N 

O 

(N 

-t— > 
o 

O: 
in 



a- 



> 

00 
00 

d 

(N 

> 

• I— I 

X: 



Yuting Ping,' El Brendon W. Lovett,^ ' Simon C. Benjamin,'-^ and Erik M. Gauger^ '-Q 

' Department of Materials, University of Oxford, Oxford OXl 3PH, United Kingdom 
^SUPA, School of Engineering and Physical Sciences, Heriot-Watt University, Edinburgh EH 14 4AS, United Kingdom 
Centre for Quantum Technologies, National University of Singapore, 3 Science Drive 2, Singapore 117543 

Coupled spin chains are promising candidates for 'wiring up' qubits in solid-state quantum computing (QC). 
In particular, two nitrogen-vacancy centers in diamond can be connected by a chain of implanted nitrogen 
impurities; when driven by a suitable global fields the chain can potentially enable quantum state transfer at 
room-temperature. However, our detailed analysis of error effects suggests that foreseeable systems may fall far 
short of the fidelities required for QC. Fortunately the chain can function in the more modest role as a mediator 
of noisy entanglement, enabling QC provided that we use subsequent purification. For instance, a chain of 5 
spins with inter-spin distances of 10 nm has finite entangling power as long as the T2 time of the spins exceeds 
270 Moreover we show that re-purposing the chain this way can remove the restriction to nearest-neighbor 
interactions, so eliminating the need for complicated dynamical decoupling sequences. 



Spin chains with nearest neighbor XY coupling mediate co- 
herent interactions between distant spin qubits with fixed lo- 
cations, and can thus serve as channels to transfer quantum 
information . An important application would be to in- 
terconnect distant sub-registers of parallel parts in a scalable, 
solid-state quantum computer |4], e.g. in a diamond-based ar- 
chitecture at room temperature USD. With an observed room- 
temperature coherence time of 1.8 ms |@], the electron spin 
of individual nitrogen-vacancy (NV^) defects in diamond is 
a promising candidate for a qubit I?]: Initialisation, coher- 
ent manipulation and measurement with nanoscale resolution 
(~ 150 nm) have already been experimentally demonstrated 
using optical techniques under ambient conditions jsfl- In ad- 
dition, the long-lived ^^N nuclear spin (/ = 1/2) associated 
with each NV^ center can act as a local, coherent memory, 
accessible via the hyperfine coupling ||^[l3|. A universal set 
of quantum operations between the nuclear memory spin and 
the processing electronic spin qubit within each NV~ center is 



available with microwave and radio-frequency pulses 111 1I4I3I1 . 
Two NV~ centers with only a small separation (r < 10 nm) 
may be entangled through direct electron spin dipole-dipole 
coupling as long as a T2 time on the order of milliseconds 
can be maintained 11411 . However, individual addressability of 
the NV^ center qubits demands larger separations of several 
tens or hundreds of nanometers and the direct interaction 
becomes too weak. 

A recent proposal (isl] suggested a chain of N implanted ni- 
trogen impurities (each with a "dark" electronic spin- 1/2) as a 
coherent quantum channel to transfer quantum states between 
distant NV~ centers at room-temperature (see Fig.[T^). Here, 
the electron spins of the NV~ centers and the nitrogen impuri- 
ties interact with each other through nearest-neighbor dipole- 
dipole coupling iflH 171. Importantly, the scheme does not 
require individual control of the chain spins, instead relying 
on global resonant driving fields to turn the effective Hamil- 
tonian into an XY exchange model 11211 . Reliable quantum 
state transfer (QST) between distant NV~ centers can then be 
achieved through this spin-chain channel ifTsll . 

An important question is how well the above described QST 
performs in the presence of realistic errors and imperfections. 
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FIG. 1. a) Two distant NV^ centers are coupled through a spin chain 
consisting of nitrogen impurities; all spins are assumed to interact 
through nearest neighbor dipole-dipole couplings. The intra-chain 
and qubit-chain coupling strengths k and g, respectively, are deter- 
mined by the geometrical arrangement of the defects, and the latter 
can be tuned through the NV^ center ground state structure to en- 
able high fidelity state transfer for g <ti k/VN fisll . b) Illustration 
of the equivalent fermionic tunnelling picture: the channel possesses 
TV energy levels with spacings \E„ — En±i\ ~ k/N. When tuned 
into resonance with a channel level E„, an NV^ spin excitation tun- 
nels through the channel and emerges on the opposite NV^ center 
after a time r„ — 7r/(\/2r„), where r„ is the effective tunnelling 
rate. Off-resonant coupling to other levels is negligible provided that 
r„ <^ \En — En.±i\ which is equivalent to g k/^/N (see text). 



Since the QST is based on free Hamiltonian time evolution, 
one would expect that timing precision and low levels of cou- 
pling disorder (due to defect implantation imprecision) are 
crucial ingredients for achieving a high fidelity llSll . How- 
ever, for the regime g ^ h/VN, QST between the NV^ cen- 
ter spin qubits through the channel has been found to be sur- 
prisingly robust against disorder in the intra-chain coupling 
K 1 19(1. The dominant error for the QST then arises due to 
the inevitable environmental decoherence of the spins in the 
channel. Focussing on such quantum noise, we argue in this 
Letter that this is indeed the limiting factor, and that current 
diamond-based architectures will fail to meet quantum error- 
correction thresholds whilst still being able to support dis- 
tributed quantum information processing. 

We consider the following effective Hamiltonian for a chain 
consisting of N spins plus two register spins located at either 
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end of the chain (see Fig.[T]l: 
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where cr!|_ = (cr^ ± i(Ty)/2 are the PauH operators act- 
ing on spin i, the coupHng strengths k and g are pro- 
portional to the cubed inverse inter-spin-distance and 
H.c. denotes the Hermitian conjugate. This Hamiltonian 
can be realised by applying global resonant microwave fields 
-ffdrive — '^x COS ujit with appropriate intensities Hi 

in the presence of a constant magnetic field applied in the z- 
direction. The Uj denote the spin energy splittings including 
the Zeeman and hyperfine components ulM (see Supplemen- 
tary Information 13811 for more details). Importantly, the basis 
of the above Hamiltonian is rotated from the physical basis 
according to {x,y,z) — ?• {z,~y,x) |[l2, 38|. Therefore, a 
T2 (Ti) process acting on the physical spin corresponds to a 
spin-flip (phase-flip) error in the basis adopted for the Hamil- 
tonian ([TJ. 



Following Ref. H15I1 we proceed by applying the Jordan- 
Wigner transformation to Eq. ([T]i, yielding 
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with fermion creation and annihilation operators 

cr'+ exp[ - i7rX;*Io '^W-]' := exp[i7rX;*Co '^W-] 
which observe the anticommutator relations {ci,Cj} = dij 

and {ci,Cj} = = {c|,c|} ll20ll . This transforms the spin 
degree of freedom into the presence or absence of a fermion 
at the relevant system site. 

For g K the coupling of the end NV~ qubits to the chan- 
nel H' can be treated perturb atively lIlSll . Diagonahsing the 
first term of Eq. (|2]i (and H.c.) yields 
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allowing us to write the second terms of Eq. (|2]i as follows 

Tl ^Cq/„ + ( 1) C^_|_j^/„ 
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where /t = /^W- 



N+l 



for n = 1,2, .. 
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and r„ = 
N liMl. The 



ff/V(A^ + l)/2sin- 
channel is now described by iJp and possesses N modes with 
energies £"„ ~ 2 k cos 7^^. One can tune the energy of both 
NV~ centers to £"„ by means of applying an appropriate de- 
tuning to 5"°'^+^ in the rotated basis. This couples the NV^ 
centers resonantly to the channel mode n with a tunnelling 
rate r,i (see Fig. [TJ)). Based on pure Hamiltonian evolution, 
the first full swap of the two end fermions occurs after a time 
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The final 'swapped' state acquires a controlled phase depend- 
ing on the total number of fermions in the system jlSi 12311 . 
Since for g <^ k/\/N off -resonant coupling to other chan- 
nel modes is negligible, the phase arises from the single mode 
and, while generally unknown at finite temperature, is well- 
defined; it can then be corrected by employing a two-round 
protocol Ill2lll5ll24ll38ll . Based on this protocol, single-mode 
coupling QST is independent of the initial chain state. 

It is readily seen that the tunnelling rate r„ reaches its max- 
imal value for odd and 71 = (A^ + 1) /2 with E,, = 0. In 
this case, no detuning is required, and the end NV^ qubits 
resonantly couple to the zero-energy mode of the channel. We 
adopt this case of odd N as the ideal implementation of the 
protocol. Since our calculations are not restricted to the single 
excitation subspace used in Ref. lfl9ll . we specifically consider 
the two cases of iV = 3 and N ~ 5 resulting in simulations 
involving 6 and 8 spins (chain plus two registers and one an- 
cilla) as explained further on. We simulate the full system 
dynamics by numerically integrating a Lindblad master equa- 
tion I25II using Hamiltonian Unless stated otherwise, we 
take the inter-spin spacing in the chain to be rN,N = 10 nm 
(k = 52 kHz) while g remains fully tuneable. 

Numerous studies seek to determine a threshold below 
which scalable quantum computing is in principle possible. 
Typically, papers quote the threshold obtained by equating 
the error rates in state preparation, measurement, and qubit- 
qubit operations (generally, the latter are the most crucial). 
For qubits embedded in diamond and linked by spin chains, 
the most relevant thresholds are those for architectures where 
qubits are arranged on a lattice-like structure with each qubit 
being 'wired' to only a few others. For this case, the threshold 
error rate is of order 1%, with 1.4% having been obtained in 
certain circumstances ll26ll . Here we will take the target error 
rate to be 0.75% (i.e., a fidelity requirement of 99.25%), asso- 
ciated with the widely studied topologically protected cluster- 
state approach 12711 . Note that to avoid diverging resource re- 
quirements, one would not wish to build a computer with per- 
formance near the threshold; practically one might target an 
error rate ten times below the threshold [28]. The fidelity of 
99.25% therefore errs on the permissive side as a requirement 
for real- world QIR We demonstrate in the Supplementary In- 
formation 13811 that this requires a physical T2 time of 25 ms 
for the N = 5 chain when using an optimally tuned ratio of 
g/K By contrast, the longest measured T2 time of the 
NV~ center is an order of magnitude below this number at 
1.8 ms 10]. Reported coherence times of the nitrogen impurity 
are much shorter (e.g. 5.5 \Jis at room-temperature and 80 [is 
at 2.5 K ^2^ ). While this should be improvable, the nitrogen 
impurity is unlikely to substantially surpass the NV^ center 
In reality, chains longer than N ^ 5 will be desirable to prop- 
erly separate the NV~ centers. We therefore conclude that 
using QST as a fault-tolerant two-qubit operation may well be 
an infeasible target for any foreseeable technology. 

Nonetheless, the suggested spin-chain quantum bus may 
still be able to support the distribution of entanglement, open- 
ing the possibility of employing distillation protocols to ere- 
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ate high-fidelity entanglement over several runs BOll . Re- 
cent studies show that highly imperfect inter-site links can 
be tolerated given that one has three or four qubits at each 
local site 131 
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The key enabling property of the chan- 
nel is then simply that it should be 'quantum' in contrast to 
a 'classical' channel, i.e. a channel that is capable of trans- 
mitting a finite amount of entanglement with each run. Iden- 
tifying the transition from quantum to classical channel with 
respect to realistic decoherence processes is the main purpose 
of this Letter. To address this question we attempt to transfer 
one half of a Bell state through the channel. More specif- 
ically, the near (i = 0) NV~ center starts off in a maxi- 
mally entangled singlet state with an additional ancilla spin, 
14--) ^ (|0)^|l)o- |l)^|0)o)/v^ and our observable is the 
entanglement of formation Ep 13311 between this ancilla and 
the remote NV~ center at position iV + 1 ll40ll . 

Let us start with TV = 3 as the shortest nontrivial odd chain 



1 4111 . Fig. lit shows that this channel is robust against physi- 
cal spin flip (i.e. Ti-type f4^) errors, which act as phase-flip 
errors in the basis of Eq. ([T]i. A high degree of transferred en- 
tanglement is achievable for Ti > 1 ms, and a finite amount 
of entanglement survives in the presence of much larger error 
rates. Longer chains also remain robust against this type of er- 
ror; for the single excitation subspace, we have simulated odd 
chains up to iV = 21 obtaining similar qualitative results 13811 . 
Further, we have confirmed that this behavior is independent 
of the initial chain state, and that our results are in agreement 
withRef. (H. 
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FIG. 2. Entanglement of formation Ef between the ancilla and the 
remote NV^ center qubit as a function of the transfer time r under 
an independent (a) spin-flip and (b) phase-flip error model for the 
nitrogen spins. The errors are applied as Lindblad operators on each 
chain spin: (a) Li = tj^. with rate 7 = l/Ti; (b) Li — o^i with 
rate 7 = I/T2. The initial state state is !*")„o |000) |0)j^_|_^, and 
g = K/{lQy/N) with k = 52 kHz. 



On the other hand, dephasing of the physical spins, charac- 
terised by their T2 time, induces spin-flip errors in the compu- 
tational basis, with a much more damaging effect on the entan- 
gling power of the channel [see Fig.|2j3]. For weak coupling 
between the NV^ centers and the channel, the QST becomes 
classical, i.e. Ep vanishes completely, for T2 < 1.6 ms. 
Applying the same model to the N = 5 chain, the tran- 
sition between a classical and a quantum channel occurs at 
T2 = 3.8 ms. This threshold increases further for larger N 
due to the transfer time becoming longer [see Eq. (|4|i], as ef- 



fectively more spins are exposed to the noise for longer 

Unlike the high fidelity transfer of a particular quantum 
state, entanglement distribution does not benefit from cou- 
pling to only a single channel mode, allowing us to explore 
stronger coupling by varying g. Fig. [3^ shows the T2 thresh- 
old which separates the quantum from a classical data bus for 
chains with TV = 3 and TV = 5, and we see that the highest 
tolerable dephasing rate occurs in the interval g e (0.9k, k). 
Expressed as a T2 time threshold this gives 120 \JlS (270 \Jis) 
for the TV = 3 (TV = 5) chain with k = 52 kHz. 
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FIG. 3. (a) T2 threshold separating classical (Ep — 0) and 
quantum (Ep > 0) state transfer. The initial chain states are 
l*")„o 1000) 10)^+, and jOOOOO) |0)^+i, respectively, but 

the threshold for other initial states is similar, (b) and (c) Maximally 
achievable entanglement of formation Ep between the ancilla spin 
and the remote NV~ center spin as a function of g/n and indepen- 
dent channel spin flip rate 7=1 /T2 for TV = 3 and TV = 5. The 
channel is entirely classical in the white regions. The inset in (b) 
shows the variation in Ep occurring for the eight different initial 
states of the channel spins in the TV = 3 case; in all cases the range 
of values of Ep is small and always less than 0.07. Random checks 
for the TV = 5 case suggest a similar behavior for longer chains. 



In Fig.[3]3 and we show the entangling capacity of the 
TV = 3 and TV = 5 channels as a function of g and the de- 
phasing rate 7 = 1/12. We note that high channel entangling 
power is only realised for a small decoherence rate of the ni- 
trogen spins, even for larger g. Unsurprisingly, the channel is 
highly inefficient close to the classical threshold. 

To 'connect' two NV^ centers separated by a fixed distance 
of 40 nm, we consider a chain comprising three and five nitro- 
gen spins. In both cases the register spins are assumed to be 
10 nm away from the ends of the chain, so that the chain's 
inter-spin distances are rN,N = 10 nm and tn^n = 5 nm, 
respectively. The first case then corresponds to the TV = 3 
chain we have studied so far, but for the TV = 5 case we 
obtain a minimal T2 ~ 34 |is for quantum communication 
(for g K, k) which is significantly shorter than 120 |j.s for the 
TV = 3 chain. The benefit of stronger coupling therefore out- 
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weighs the drawback of a longer chain for a fixed separation 
between the register spins by speeding up the transfer process. 
However, how many impurity spins can be deployed will cru- 
cially depend on the achievable implantation precision, since 
even sub-nanometer imperfections entail significant coupling 
disorder for closely spaced chains. 

So far, we have assumed nearest-neighbor couplings [see 
Eq. (H iL which could be realised through dynamical decou- 
pling il2ll . In practice, one needs a coupling strength 'high 
pass filter', which would typically reduce effective coupling 
strengths and add further complexity to experimental imple- 
mentations. We now relax this assumption and include all 
pairwise couplings with a dipole-dipole interaction \/r^ dis- 
tance dependence, including register to chain spin couplings. 
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FIG. 4. Top: Ef between the ancilla spin and the remote NV~ cen- 
ter as a function of the transfer time r for all eight initial chain states 
of the A'^ = 3 non-nearest neighbor chain. All spins are assumed 
equally spaced, with the coupling strength for nearest neighbors be- 
ing 52 kHz. Bottom: Ef for the initial state |010) 
with decoherence processes of independent phase-flip noise (left) 
and spin-flip noise (right) in the basis of Hamiltonian ([T} correspond- 
ing to physical Ti and T2 processes, respectively. 



Fig. m demonstrates that the entanglement transfer for 
a uniformly spaced spin chain now depends on the initial 
chain state. However, the first maximum of all initial states 
coincides and gives rise to a reasonable degree of trans- 
ferred entanglement. Thus we pick the initial chain state 
I^~)a0 I^IO) |0)jY+i as an example that is indicative of the 
performance to be expected. Interestingly, Ti and T2 type 
processes degrade the Ep with much more similar severity 
compared to the nearest neighbor only coupling case 14311 . 
Nonetheless, there still exists a threshold in the T2 coherence 
time, and to ensure quantum communication for an arbitrary 
initial state, we require T2 > 140 ixs, which is only slightly 
longer than for the nearest neighbor chain. For the case of the 
N = 5 chain (not shown), we have determined the transition 
threshold to be T2 > 340 [xs, which again only represents a 
modest increase. Taking into account the complexity and ex- 
pected reduction in effective coupling strengths from dynam- 
ical decoupling sequences, it may thus be more practical to 



drop the restriction of only nearest neighbor interactions. As 
we show in the Supplementary Information, the state transfer 
still remains similarly robust to small coupling-strength disor- 
der il. 



In conclusion, by employing numerical simulations we 
have studied the impact of inevitable decoherence processes 
on the entanglement capacity of a spin-chain bus that is re- 
alised by dipole-dipole coupling of crystal defects. Limiting 
our discussion to chains of length iV = 3 and TV = 5 as likely 
candidates for a first experimental demonstration (and also for 
numerical tractability) has allowed us to obtain insight into 
important characteristics of such a protocol in the presence 
of realistic noise. Our conclusions will be equally relevant 
for longer chains, which are known to be even more suscep- 
tible to decoherence 134143711 . We have shown that directly 
meeting quantum error correction thresholds remains infeasi- 
ble even if all spins possessed the exceptional coherence time 
of NV^ centers. In contrast, the distribution of a finite amount 
of entanglement appears realistic with current systems, offer- 
ing the possibility of applying distillation protocols to boost 
the transmitted entanglement with additional local operations 
and classical communication jsoll . 
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In this Supporting Information document, we present further information and calculations supporting the 
conclusions of the main Letter. In particular, we discuss the level structure of the NV^ center, and give an 
explicit derivation of the XY spin chain Hamiltonian used in the main text. Furthermore, we show that fault- 
tolerant two-qubit gates require unrealistically long T2 times. Finally, we present our results for Ti processes in 
longer chains and the effects of coupling strength disorder. 



ELECTONIC SPIN QUBIT OF NV DEFECTS 

The nitrogen-vacancy colour defect in diamond consists of 
a substitutional nitrogen atom and an adjacent vacancy. In its 
negatively charged state, the NV~ center traps an excess elec- 
tron and possesses a paramagnetic ground state {S = 1) with 
extraordinarily long spin lifetime. For each NV^ center, the 
spin-triplet ground state '^A consists of the nis = and the 
degenerate (by Csi, symmetry) rris = ±1 sublevels, split by 
2.88 GHz |llj,|2|]. An external magnetic field lifts the degener- 
acy between the rris = +1 and the nis = —1 sublevels. The 
637 nm optical transitions between the ground and the excited 
■^E triplet states are predominantly spin-conserving 0], 
but an intersystem crossing (ISC) occurs via the singlet ^A 
state Hill (see Fig. [IB- 




ms = 



FIG. SI. Simplified electronic structure of the NV^ colour center. 
The zero-field splitting between the rris = and the rUs = ±1 
ground triplet is 2.88 GHz, and transitions between different spin 
levels can be effected through resonant microwave pulses [UsQ]- 
The singlet level is metastable and provides a route for intersys- 
tem crossing relaxation (red arrows) from the excited triplet to the 
ground state: Thick red lines have associated ISC rates that are sev- 
eral orders of magnitude faster than the thinner red line, and dashed 
transitions can be considered to be negligible (HQ]. 



EFFECTIVE SYSTEM HAMILTONIAN 

Following Ref. ^ we derive the effective system Hamilto- 
nian for the chain shown in Fig. 1 of the main text. We assume 
the presence of a constant magnetic field of strength B in the 
z-direction which is defined by the symmetry axis of the NV~ 
center, i.e. the [111] crystal axis. The Hamiltonians for indi- 
vidual NV^ and N defects are given by |[l3.llill, respectively, 



+i?(5r)'+^Nvi-s'^^ 

^^N = gefJ-sBS^ - gnUnBIz + S • An 



I, 



(SI) 
(S2) 



where ((?«) is the electron (nuclear) g-factor, iib if^n) is 
the Bohr (nuclear) magneton, ^nv (An) denotes the hyperfine 
coupling constant (tensor), D is the zero-field splitting for the 
NV^ center, and S^'^^) (I) is the full electronic (nuclear) spin 
operator. 

For each NV^ center, we encode the qubit basis states |0) 
and |1) in the rris ~ and rris = +1 sublevels of the ground 
triplet (see Fig. lSlI ), respectively. Expressed in the computa- 
tional basis of the qubit, we can approximate Eq. ( ISlb as: 



rrqubit 

ll\]\f — 



A. 



(S3) 



where the electronic Zeeman energy :^ D + gfjigB, 
and az is the usual PauU z-operator acting on the qubit. The 
effective Hamiltonian above does not include the heavily sup- 
pressed hyperfine spin-flip terms, nor the negligible nuclear 
Zeeman term (^nv, .gnMn-S ^ ^0 ^) 

For the z*'* nitrogen defect with S = 1/2, Eq. ( IS2b can be 
similarly approximated as 



To use the electron spin of an NV~ center as a quantum 
bit, one can, for example, encode the states |0) and |1) into 
the rris = and nig = +1 sublevels of the ground triplet, 
respectively. Optical pumping then polarizes the qubit into the 
|0) state, and resonant microwave pulses enable single qubit 
operations, such as the creation of a coherent superposition 
state (|0) + |l))/\/2 iH 01. The spin can also be read out 
optically with a fluorescence technique 10, • 



-Hn = yCTz + — J^CTz, (S4) 

where ujq := g^fisB w 10 GHz, and al acts on the spin qubit 
of the i*^ N impurity. The hyperfine couplings AL| depend on 
the Jahn-Teller orientation of each N defect 191 1 1 2ll 1 311 . 

The magnetic dipole-dipole interaction between two elec- 
tron spins i and j (both N defects, or one N and one NV~ 
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center) is generically given by 



T2 THRESHOLDS FOR FAULT TOLERANCE 
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(S5) 



where /io is the vacuum magnetic permeability, r is the sep- 
aration between the spins, and f denotes the unit vector con- 
necting the two spins, here assumed to be parallel to the x- 
direction (as is the case in Fig. 1 of the main text). In the last 
step, we have neglected the spin-flip terms as before, since 
A'o.9eMB/(^''^) — 52 kHz < ^nv,-AL|| for a spacing of 
r = 10 nm Jglll^l. We can therefore 13 write the Hamilto- 
nian for a pair of spins as 



(wo + 5i) 



(S6) 



i=l,2 

Hn.n\ ^ g (Tz(^z^ 7, (^z^ n CTz' (^') 



where the hyperfine terms 5 := ±Anv/2 and Si := ±^^||/2 
(i = 1,2, N); the dipolar coupling strengths are k := 

Mo5eMl/(7r''N,N) and 5 := Mo5?Ml/(7r''N,Nv)- 



Eqns. ( IS6I 1 and ( IS7b are easily generalised to a nearest 
neighbour coupled chain. We assume that the entire chain is 
driven by the following resonant global fields 



N 



H drive = ^ Cr* COS UJ it 

i=l 

+ ^ r^o cos(c< 

j=0,N+l 



6)t, 



(S8) 



We simulate the quantum state transfer (QST) for spin 
chains with odd N in the weak coupling regime, g = 
k/{10\/N), and for a 10 nm inter-chain spacing. After the 
transfer the acquired phase ±1 of the target state is corrected 
to enable a direct comparison with the initial state (in practice 
the phase would be cancelled by employing a two-round pro- 
tocol 10 15, 3l). To evaluate the fidelity of the transfer pro- 
cess, we use the measure F^{p, a) = Tr {y/ ^/p(Jy/p) ifisll . 
The N defect spins are subjected to independent physical T2 
processes, realised as spin flips in the computational basis. 

Fig. [S2] shows that meeting the well-studied topological 
fault-tolerance threshold > 99.25% ()3] requires unre- 
alistically long coherence times of the defect spins. More 
specifically, the shortest non-trivial TV = 3 chain achieves a 
sufficiently high fidelity only for T2 > 70 ms whereas this 
number increases to T2 > 350 ms for iV = 5. 
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FIG. S2. Fidelity F of the transferred state p with respect to the in- 
put CT = 1+) (+|, through the A'^ = 3 chain as a function of the trans- 
fer time T. The initial state is \+)q |000) |0)^^]^, and we operate in 

the weak coupling regime g — K/(10A/iV) for k = 52 kHz. The ni- 
trogen spins experience independent dephasing at a rate 7=1 /T2 . 
Whilst the curves seem to coincide in the main plot, the inset shows 
a zoomed-in view near the maximum on the scale relevant for fault 
tolerance. 



where all four possible frequencies uji ~ ujQ + Si are applied to 
address the nitrogen impurities 1^. The field intensities flo, ^ 
are chosen to fit in the hierarchy k ^ 51, fio ^ Wi, + 6. 

Making a rotating wave approximation in the usual ro- 
tating frame (see Appendix) and adopting the rotated basis 
(x, y, z) — > {z, —y, x) |@] then yields an effective XY model 
for the chain: 



N-l 
i=l 

+ ^ g{alaL+'+aW++'), (S9) 

j=O.N 



where ct^ = « ± zct;)/2. 

Eq. ( IS9l l is the effective system Hamiltonian given as 
Eq. (1) of the main text. Controlling the magnitude flo can 
also effectively tune the coupling q between the NV^ center 
spins and the N defect spin chain 



As described in the main text, in the weak coupling regime 
NV^ center excitations tunnel through the (single) zero- 
energy mode of the chain, and off -resonant coupling to other 
modes is negligible. This enables high-fidelity quantum state 
transfer for long enough coherence times. However, because 
transfer time t ~ 1/g is longer for weaker g, the state transfer 
is more susceptible to decoherence, and it may thus be ad- 
vantageous overall to use stronger g at the cost of coupling to 
several modes and sacrificing some theoretical fidelity. In the 
following table we list the T2 times in milliseconds required 
for achieving an error rate below 0.75% for different g/n ra- 
tios with K = 52 kHz: 



5/« 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


1 


N = 3 


41 






28 








4.5 


7.2 




N = 5 


160 


110 


83 








91 


25 







where '-' denotes that the desired fidelity is never achieved 
for the initial state |+)o |000) |0)^^j^. The fluctuating behav- 
ior seen in this table is consistent with the results reported 
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in Ref. d. We conclude that for the studied chains fault- 
tolerant quantum computation demands coherence times of 
several milliseconds for the = 3 chain and a few tens of 
milliseconds for the more interesting = 5 chain, even when 
dropping the weak coupling constraint and in the absence of 
any other imperfections. 



Ti PROCESS IN LONGER CHAINS 

In the main text, we show the effects of Ti and T2 pro- 
cesses for chains of length = 3 and 5. Considering only Ti 
processes and restricting ourselves to the zero and single exci- 
tation (computational) subspace ll20ll significantly reduces the 
numerical complexity, allowing us study longer chains. 
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FIG. S3. Entanglement of formation Ef between the ancilla spin and 
the remote NV~ qubit as a function of the transfer time r under an 
independent, physical bit-flip model for nitrogen spins with rate 7 = 
l/Ti. The chain length is A'^ — 15, and all chain spins are initially 
in the state 1 0) . Only nearest neighbour interactions are included. 



Fig.lS3lillustrates that a sizeable amount entanglement can 
be transmitted through a A^ = 15 chain for Ti times as short 
as a millisecond, both for strong and weak coupling of the 
NV^ centers to the chain. When the bit flip rate gets small, 
Ti > 0.1 s, the weak coupling approach possesses a much 
higher entangling power, despite taking substantially longer. 
For Ti 10 ms fl9ll . however, both approaches are com- 
parable. Taking into account the increased robustness of the 
strong coupling case against T2 processes (see the main text) 
suggests that overall 5 w k is likely the better choice for tack- 
ling decoherence. 

Fig. IS4I considers different chain lengths with g = k cou- 
pling, showing the optimal transfer duration t and maximally 
achievable Ep for each case under an independent bit flip 
model with Ti = 10 ms. Unsurprisingly, the performance 
of the chain is more affected for longer chains, although the 
reduction is much less drastic than for T2 noise, and a fi- 
nite amount of entanglement can be transferred even for long 
chains. 



COUPLING-STRENGTH DISORDER 

In this section, we simulate disorder in the inter-chain cou- 
pling strength k, as would arise from imprecisions when im- 
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FIG. S4. Plot of the transfer time r (the first peaks in plots like 
Fig. lSBT a)) and the resulted Ef values between the distant spin and 
the (7V+1)*''NV" qubit as a function of the channel length TV, with 
fixed inter-chain spacing rm = lOnm {g = n — 52kHz), under an 
independent phase-flip model (in the rotated basis x -tr^ z) on each 
nitrogen spin with practically relevant rate 7 — 1/Ti — lOOHz 
(Ti = 10ms). The channel spins are initially all in the state |0), and 
interact with their nearest neighbours only. 



planting the nitrogen impurity. For numerical convenience, 
we once more restrict our calculations to the single excitation 
subspace. We assume the spacings between the neighbour- 
ing spins obey a Gaussian distribution around the mean value 
tnn = 10 nm (k = 52 kHz). Our results are averages over a 
hundred independent runs for each data point. 
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FIG. S5. Maximally achievable Ef between the ancilla and the 
remote NV^ qubit for different chain lengths N. The inter-chain 
spacings are assumed to follow a Gaussian distribution with mean 
^NN ~ 10 nm with a standard deviation corresponding to 5% dis- 
order (~15% disorder in the inter-chain coupling strength k). The 
channel spins are initialised in the state jO), and each data point is 
the result of an average of 100 independent runs. The error bars in- 
dicate a 95% confidence interval. The 'nearest neighour' case uses 
the more practically relevant strong coupling regime g — Kmean- For 
'beyond nearest neighbour coupling', all pairwise spin couplings are 
included in accordance with the dipolar 1 /r^ distance dependence. 



Fig.lSSIshows that even a sizeable 15% spread in the k dis- 
tribution does not have a catastrophic impact on the entangling 
capacity of the chains. In fact, the reduction of achievable Ep 
is rather similar to that obtained from a Ti time around 10 ms. 
While the difference is not huge, the nearest-neighbour cou- 
pled chain proves consistently more robust. 

However, for the short chains considered in the main text. 
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a small amount of disorder in k is unlikely to be the limiting 
factor preventing entanglement distribution. Ultimately, we 
expect the major challenge for this protocol is attaining suffi- 
ciently long T2 times to overcome the limitations discussed in 
the main Letter. 



nian ( IA3l l then becomes 



12 ^ 



1=1,2 
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APPENDIX: ROTATING WAVE APPROXIMATIONS 



Under the additional driving field of Eq. ( IS8b . the total 
Hamiltonian for two nitrogen electronic spin qubits reads 



1=1,2 
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The unitary transformation for moving into the rotating frame 
is given by 
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where 0i 



(4) 



and 6' 



2 — " -2(3) - V y 2 

ing the transformation Hrp = iUW + UHU^ and making 
the usual rotating wave approximation (RWA), justified since 
il,K uji, yields: 



/ 



(A2) 



Apply- 



RWA 



'4 ~K f 

f f 



j=l,2 



In the rotated basis with (a;, y, z) 



(A3) 



i.e. |0(1)) 



-(-)> = (|0) ± |1))/V2), the Hamilto- 



:^ . + ala\) + ^a^ (A4) 

•i=l,2 

Here, a second RWA was made in the last line by neglecting 
the non-spin-conserving terms, which is valid when k <C 

The same procedure can be generalised to a longer chain 
straightforwardly, allowing us to arrive at the desired nearest- 
neighbour interaction Hamiltonian i 
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